Low-energy subgap states and the magnetic flux periodicity 
in d-wave superconducting rings 
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Wave functions of low-energy quasiparticle subgap states in d-wave superconducting rings, 
threaded by Aharonov-Bohm magnetic fiux, are found analytically. The respective energies are 
closest to the midgap position at small magnetic fluxes and deviate from the Fermi surface due to 
the Doppler shift, produced by the supercurrent. The Doppler-shifted zero-energy states result in a 
paramagnetic response of the ring at small fluxes. The states exist only for even angular momenta 
of the center of mass of Cooper pairs, in agreement with recent numerical studies of the problem. 
This macroscopic quantum effect in d-wave rings results in broken h/2e periodicity, retaining only 
the h/e periodic behavior of the supercurrent with varying magnetic flux. 
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The h/2e periodic dependence of the current on the 
magnetic flux in superconducting rings and hollow cylin- 
ders is considered frequently as an inevitable consequence 
of the electronic pairing in superconductors. This point 
of view is fully supported by the Ginzburg-Landau the- 
ory, which contains only the Cooper pair charge 2e and 
always predicts the magnetic flux period h/2e. Numer- 
ous experimental results seem to be also in favor of such 
conclusions, including first observations of the h/2e flux 
quantization in hollow cylinders P, Q , the Little-Parks 
effect and the flux quantization of Abrikosov vortices 
0, 0|- However, the Ginzburg-Landau approach applies 
on the scale much greater than the Cooper pair size and, 
strictly speaking, at temperatures near Tc. It is, in gen- 
eral, not applicable to mesoscopic rings at low temper- 
atures, where a microscopic approach, which is not so 
plain with respect to the problem 0, 0, ll], has to be 
explored. 

According to the microscopic BCS theory, normal- 
metal electrons with projections hM and hAI of their 
orbital angular momenta along the ring axis, form in s- 
wave superconductors Cooper pairs with the angular mo- 
mentum hq = h{M + M) of their center of mass 0, [^. 
Respectively, a Bogoliubov quasiparticle in the super- 
conducting ring represents a superposition of an electron 
with the angular momentum hM and a hole with —hM, 
and the difference is hq. Since M and M have to be 
integers and can be represented under the conditions in 
question as M = {£ + q)/2 and M = {-i + q)/2, the 
quantities £ and q can take simultaneously either even 
or odd integral values [1, il. The key point is that the 
single valuedness of the wave functions is protected by 
the quantities {£ ± q)/2, whereas only q/2 enters the ex- 
pression for the supercurrent via a standard combination 
g/2 — 4>/<&o with the magnetic flux. As a result, the fac- 
tor g — 2<I>/$o arises in the microscopic BCS theory and 
leads to the magnetic flux quantum $o/2 in superconduc- 
tors^a half of the normal-state flux quantum $0 = h/e 



The condition of the single valuedness of the wave 
function, which justifies the presence of the supercon- 
ducting flux quantum <&o/2, does not necessarily result 
itself in the h/2e periodicity, however. An important 
additional condition, which has to be satisfied for en- 
suring the periodicity, is the degeneracy of the respec- 
tive states. In general, the degeneracy takes place be- 
tween states of a superconducting ring pierced by mag- 
netic fluxes, which differ by an integral number of $0i i- 
e. by an even number of the superconducting magnetic 
flux quanta ^o/2 = h/2e. A difference by odd numbers 
of the quanta is, in principle, physically distinguished 

0| • This circumstance represents a signiflcant interest 
since it could lead to a ft,/e-periodic component in behav- 
iors of superconducting rings, which removes the specific 
superconducting h/2e periodicity or, at least, makes it 
approximate. 

On the other hand, there are microscopic arguments, 
which substantially restrict possible deviations from the 
h/2e periodicity and explain from the microscopic point 
of view its numerous experimental observations in super- 
conducting rings. The arguments give also a general idea 
for further search for the conditions, when the breaking 
of the periodicity is noticeable. As this follows from the 
BCS theory, the breaking could become observable, if, in 
carrying out the statistical averaging, one cannot always 
replace the sum over states near the Fermi surface by the 
respective integral. For this reason the discreteness of the 
angular momentum is of crucial importance for making 
it possible to distinguish between the states of a ring, 
which differ by odd superconducting fiux quanta. For s- 
wave pairing, this can happen in mesoscopic rings with 
characteristic sizes of the order of the coherence length 
^0 or less Among various possible complications of 
the standard BCS approach which arise in rings of such a 
small size, probably, the main one is that the discreteness 
of the states and the pair breaking effect of the supercur- 
rent can destroy the superconductivity there, at least, 
within some range of the magnetic fiux fl^ . This com- 
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plicates an experimental observation of the magnetic flux 
periodicity in a superconducting state of the mesoscopic 
s-wave rings. The fluctuations of the type of quantum 
phase shps, as well as effects of the Coulomb blockade, 
could also play an important role there. 

The problem has been investigated recently in RcfH 
for d-wave rings at low temperatures, where a noticeable 
violation of the h/2e periodicity, associated with discrete 
current-carrying low-energy states, was identified. One 
of the most striking results of Ref@, obtained within 
numerical selfconsistent calculations, is that in d-wavc 
superconducting rings quasiparticle subgap states with 
quite low energies exist only for even q. By contrast, a 
considerable empty spectral gap appears in d-wave rings 
with odd q. In the presence of low-energy states, a com- 
paratively small variation of the Dopplcr shift with the 
magnetic flux turns out to be sufficient for generating an 
alternation of the diamagnetic and paramagnetic states 
in the ground state of the d-wave ring. As a result, the 
breaking of the h/2e periodic dependence of the supercur- 
rent on the magnetic flux has been found to be observable 
in d-wave rings, whose size can significantly exceed the 
coherence length ^o- 

In the present paper an analytical approach will be de- 
veloped for studying the effect, based on a comparatively 
simple model of the d-wave superconducting ring. Ener- 
gies and wave functions of quasiparticle current-carrying 
states, lying closest to the Fermi level at small mag- 
netic fluxes through the rings, will be obtained by solv- 
ing the respective quasiclassical equations. The states 
turn out to deviate from the zero energy (the midgap 
position) only due to the Doppler shift, produced by the 
supercurrent. They form a paramagnetic current, which 
dominates the response of the ring at low fluxes. Since 
the anisotropic d-wave pairing breaks the conservation 
of the angular momentum along the ring axis, the angu- 
lar dependence of the respective quasiparticle wave func- 
tions does not reduce to a simple exponential dependence 
exp[i(£±g)(y9/2], which they would possess in isotropic s- 
wave or normal-metal rings. The probability density of 
the states has its maxima at the nodes of the order pa- 
rameter. 

The phase factors of the quasiparticle wave functions 
are also modified in d-wave rings. This results in a topo- 
logical reason for the low-energy quasiparticle states to 
survive only for even q. It will be shown, that the quan- 
tity £/2 in the phase factor is replaced by the value Mp of 
the angular momentum along the ring axis of the normal- 
state quasiparticlcs at the Fermi energy, for a given trans- 
verse channel. Since 2Mp is always even, the exponential 
exp[i(2Mi^±(7)(^/2] is compatible with the single valued- 
ness of the superconductor wave-function only for even 
q and, hence, for even £. Therefore, the Doppler-shifted 
zero-energy states exist only in even g-sectors of the pair- 
ing. This signifies that the violation of the h/2e period- 
icity in the subgap spectrum of d-wave superconducting 



rings is a macroscopic quantum effect associated with 
the discreteness of the states and the specific structure 
of their wave functions. The discreteness is of the uni- 
versal character and should be observable at sufficiently 
low temperatures, when pairings with different q in su- 
perconducting rings can be distinguished. 

Subgap states Consider a narrow strongly type II super- 
conducting ring with a negligible supercurrent-induced 
magnetic flux, where the arm width L is much less than 
the penetration depth L <^ and considerably ex- 
ceeds the coherence length. The nonlocal operator for an 
anisotropic order parameter in the ring can be described 
quasiclassically and reduced, in the main approximation, 
to a local angular dependent quantity. 

The dependence of the Bogoliubov amplitudes on the 
polar angle in continuous circular (cylindrical) rings can 
be taken in the form 



u{(p) = u{ip)cxp -{i + q)ip 
r 2 

v{ip) = v{ip) exp ^-(i-q)^ 
^ A{(p) = A{lp) cxp {iq(p) . 



(1) 



Quantities u and v do not depend on the polar angle if 
in isotropic normal metals and s-wave superconductors. 
Since in d-wave rings the angular dependence of ampli- 
tudes u{ip) and v{ip) is induced entirely by an anisotropy 
of the superconducting order parameter A(ip) and dis- 
appears in the normal-metal state, the amplitudes vary 
comparatively slowly with changing ip. Taking this into 
account, the following equations of the Andreev type for 
-11(1^9) and v{(p) can be derived: 

hvp du ( \_ , \~ 

I — TT^- h I £ — 4 — mvp^ipVs ju — A[ip)v = , 



e + £,-mvp^^Vs)i + A*{ip)u = Q. (2) 



R dip 

.hvpip dv 

I ^ — 

R dip 



Here v^ is the supercurrent velocity, which takes the 
form 2mRvg = /imin^ (g — 2$/$o); ^ is the external 
Aharonov-Bohm magnetic flux and the associated vector- 
potential is Aip = $/27ri?. The normal-state quasiparti- 
cle excitation energy ^ is taken with respect to the Fermi 
level. For a given transverse channel, ^ is a discrete 
quantity, which depends on £/2 in Eqs.(l2]), and the ef- 
fective Fermi velocity of a circular motion is defined as 
VF.ip = c^C/t^Pv|pj^' = M/2R. 

As the radius of the ring is much greater than the 
atomic scale kpR ^ 1, the quasiclassical approach ap- 
plies to a circular quasiparticle motion and, in partic- 
ular, admits the condition |vF,ip| ^ |vs|, i. e. \l\ ^ 
|(7 — 2<I>/<I>o|. For the maximum value of the super- 
fluid velocity \rrrvp^^Ys\ ~ max [|A((^)|], one gets |vs| ^ 
(^F^o)"^ V'F,ip\ ^ |vf,(,3|. Hence, linear terms in Vg, as 
well as all other terms in Eqs.([2), can be considered 
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within the quasiclassical description as small quantities 
of the first order as compared with £pu{(p). 

In s-wave superconductors, the amplitudes u, v do not 
depend on Lp and the condition for nontrivial solutions of 
Eqs. results in the quasiparticle energies 



£ = mvp^^vs ± VeTW ■ (3) 



The ground-state value of q is known to possess the 
h/2e period, changing with flux as q = int[2(f>/<l'o + 
sgn(<i>)/2], up to small finite-size corrections. Such a 
dependence is assumed observable below, for example, 
if the system is field cooled through Tc for each field 
value. In general, states with different q describe differ- 
ent metastable phases of the superconductor 0, which 
usually possess large lifetimes, even if L <C Xl and the 
magnetic flux is not really trapped. If pair breaking ef- 
fects of the supercurrent are small, the energies ([3]) de- 
pend linearly on g — 2$/$o and manifest the same h/2e 
periodic dependence. 

In the absence of the supercurrent, energies ([3]) coin- 
cide with those in homogeneous superconductors. The 
Doppler shift, described by the first term in Eq.Q, 
causes some of the states to move in the subgap region 
with varying magnetic fiux. As [q — (2<i>)(<i>o)] < 1/2 for 
the ground-state behavior of q, the energy ^ can cross 
the Fermi surface at some flux value only in small rings, 
when 4i? < hvF,<p/\^\ and the pair breaking supercurrent 
strongly modifies or even destroys the superconducting 
state in the ring. In s-wave rings of larger size, ener- 
gies of the subgap states lie comparatively far from the 
midgap position. 

Andreev equations ([2]) can be transformed to Riccati 
equation [3], which has periodic coefficients in a dou- 
bly connected d-wave sample. Numerous solutions of 
Eqs. ^ will not be discussed analytically here, as this is 
mostly the problem for numerical studies. The selfconsis- 
tent numerical solution of the respective Bogoliubov-de 
Gennes equations for the tight-binding model of the d- 
wave square loops is represented in RefQ. There are, 
however, two remarkably simple degenerate continuous 
solutions of Eqs. (21) for an intrinsically real order param- 
eter A {(p), which changes its sign and vanishes after aver- 
aging over the polar angle. These solutions are of special 
interest and can be described analytically for any partic- 
ular angular dependence A{ip). For a simple model of 
the c?2,2_y2-wave order parameter, one can write in the 
main quasiclassical approximation A (ip) = Ad cos 2(p af- 
ter applying A to ([T]). For spatially constant A^, the 
"electron-like" and the "hole-like" solutions take the fol- 
lowing form (in a more general case one should make the 
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FIG. 1: The probability density as a function on the polar 
angle for the Doppler-shifted zero-energy state in the d-wave 
ring with R = 10hYF,ip/\Ad\. 

substitution in Eqs. (g]), © A^ sin(2(^)/2 J A{(p)d(p): 




(4) 



(5) 

where /3 = R^/hvp^^p, e = mYp^^Ys and constant ampli- 
tudes Ci,2 result from the normalization. 

The energy e = zLm \vf^^Vs\ deviates from the midgap 
position due to the Doppler shift and its sign depends on 
relative directions of the quasiparticle and the supercur- 
rent circulations. If the ring size noticeably exceeds ^o, 
then |e| ^ \^d\- The origin of the Doppler-shifted zero- 
energy states is associated with the change of sign of the 
d-wave order parameter. Fig. [T] displays the probability 
density of the states, which possesses the characteristic 
fourfold structure and reaches its maximum exactly at 
the nodes ip = ±j, As for superpositions of (j4|) 

and ([5]) the fourfold symmetry reduces to two twofold 
structures, the degeneracy of the states Q and (l5])could 
be lifted within higher-order approximations. 

The applicability of the quasiclassical approach to the 
solutions (HI and ([5]) can be explicitly justified. The 
terms containing the derivatives of the first order in Eqs. 

can be estimated as {hvp^^/R)du{(p)/dip ^ Af)u{ip), 
or ^ ^{(((/j), and identified as small quantities of the first 
order. At the same time, the terms which have been 
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neglected within the accuracy of the Eqs.([2]), contain an 
additional small factor Aq/ef or ^/ej? as compared to the 
first-order terms. This concerns, in particular, the terms 
with the second-order derivatives, which can be, there- 
fore, classified as small quantities of the second-order. A 
variation of the coefhcient Hvp^ip/R in front of the first- 
order derivatives in Eqs.([2]), associated with the differ- 
ence between v^ and the respective polar component of 
the Fermi velocity vf,i^, also results in small terms of 
the second order, which are beyond the accuracy of the 
equations and have to be neglected. This circumstance 
is important for further considerations. 

Broken h/2e periodicity Wave functions of the orbital 
motion have to be single- valued, if only single-valued 
gauge transformations are used and spin is not intimately 
involved in the problem (isl. [lij . This is the case, in 
particular, for the Bogoliubov amplitudes H]), (|4]) and 
([5|). The exponentials in Eq.(IT]) are always single- valued, 
since i and q take even or odd values only simultaneously. 
Hence, the amplitudes u{(p) and ^(i^) have to be single- 
valued themselves and their phases have to acquire an in- 
teger number of 2tt after going around the loop. This sig- 
nifies that the solutions Q , ([5|) really exist only if /? is an 
integer. For the last condition to be satisfied, the quanti- 
zation of the angular momentum along the ring axis turns 
out to be of crucial importance. Indeed, close to the 
Fermi surface \{i/2) — Mp] <C \Mp\ and the quantized 
normal-state excitation energy £,{(/2) takes the follow- 
ing form: = vfAP'p " Pf,v) = ^ (^/2 " Mf)- 
Therefore, /3(€/2) = {£/2 - Mp). As Mp is always an 
integer, (3 is an integer- valued parameter for even (. and, 
hence, for even q. In the case of odd £ and q the param- 
eter (3 takes half-integer values, however. For this reason 
the Bogoliubov amplitudes H]), (j4|) and ([5]) arc not single- 
valued for odd 5, changing their sign after going around 
the loop. Hence, the Dopplcr-shiftcd zero-energy states 
(|4]) , ([5]) exist in the d-wave superconducting rings in even 
g-sectors of the pairing, whereas for odd q the states do 
not arise in the gap. 

As a rule, <^ can be effectively excluded from a theo- 
retical analysis of observable superconductor properties. 
For example, the current and other obscrvablcs are ex- 
pressed via quasiclassical Green functions, which are ^- 
intcgrated. As this follows from the present paper, the 
situation can change substantially, if ^ takes discrete val- 
ues due to the angular momentum quantization. 

Paramagnetic response The total current J induced by 
the Aharonov-Bohm field ^(^(p) = ^/2'Kp in a cylin- 
drical ring, is obtained from the relation ((5£($))^ = 
-(1/c) j^{p)5A^{p)dV = -(l/c)J(5$, where £ is the 
energy and S the entropy. Thus, in even g-sectors 
of the pairing, a contribution to the subgap current 
from the two states with ±£/2 takes the form J = 



abrupt changes associated with the Sq = ±1 transitions 
have been disregarded. Within the model, the character- 
istic value of the current eYF,ip/2'!rR, carried by the low- 
energy states, coincides with that for the normal-state 
persistent current [l7j . However, its sign is determined 
by the quantity q—^, which is specific for the supercon- 
ducting state. As a result, it is a paramagnetic current. 
The current dominates the magnetic response of the ring 
in vicinities of centers of even g-sectors, including small 
fluxes at q = 0, where energies of the states (|4]), ([5]) are 
most close to the Fermi surface. The paramagnetic re- 
sponse of the zero-energy states is known to take place 
also near surfaces of d-wave superconductors [3, [l^ . 

Conclusion Energies and wave functions of the low- 
energy quasiparticle states have been obtained analyt- 
ically in d-wave superconducting rings threaded by an 
Aharonov-Bohm magnetic flux. The states turn out to 
deviate from the zero energy due to the supercurrent- 
induced Doppler shift. They form a paramagnetic re- 
sponse of the rings at small fluxes. The Doppler-shiftcd 
zero-energy states arc found to exist only for pairings 
with even angular momenta of the center of mass of 
Cooper pairs. This macroscopic quantum effect breaks 
h/2e periodic behavior of the supercurrcnt in the ring, 
in agreement with the results of Ref@. The analytical 
approach developed in the present work, demonstrates 
explicitly that the quantization of the orbital angular 
momentum and the condition for the single valuedness 
of the superconductor wave function play the key role in 
a formation of the difference between quasiparticle sub- 
gap spectra at even and odd q. 
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the derivative of the energy over the magnetic flux. 
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